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Abstract 

We use asymptotic analysis and numerical simulations to study peak-type singular solutions 
of the supercritical biharmonic NLS. These solutions have a quartic-root blowup rate, and 
collapse with a quasi self-similar universal profile, which is a zero-Hamiltonian solution of a 
fourth-order nonlinear eigenvalue problem. 

1 Introduction 

The focusing nonlinear Schrodinger equation (NLS) 

^,x) + AV + |Vf ff V< = 0, ^(0,x) = Vo(x) € H\R d ), (1) 

where x = (xi, . . . , x^) £ M. d and A = X^=i d 2 is the Laplacian, has been the subject of intense 
study, due to its role in various areas of physics, such as nonlinear optics and Bose-Einstein Con- 
densates (BEC). It is well-known that the NLS (T1J possesses solutions that become singular in a 
finite time [SS9^] . 

In recent years, there has been a growing interest in extending NLS theory to the focusing 
biharmonic nonlinear Schrodinger equation (BNLS) 

^(t,x)-A 2 V + M 2ff V> = 0, V(0,x) =Vo(x) £H 2 (R d ), (2) 

where A 2 is the biharmonic operator. The BNLS ([2]) is called "l?-critical\ or simply "critical" 
if ad = 4. In this case, equation ([2]) can be rewritten as 

^(t,x)-AV + |V'| 8/d ^ = 0, V(0,x) =Vo(x) eH 2 (R d ). (3) 

Correspondingly, the BNLS with < ad < 4 is called subcritical, and the BNLS with ad > 4 is 
called supercritical. This is analogous to the NLS, where the critical case is ad = 2. 
BNLS solutions preserve the power (L 2 norm) 

P(t) = p(o), p = M\1 



i 
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and the Hamiltonian 



H(t) = H(0) 



H = \\Ai;\\l 



1 



2ct+2 
2o-+2- 



a + 1 



In [BAKSOO] . Ben-Artzi, Koch and Saut proved that when a is in the H 2 -subcritical regime 



jo<a d<4, 
\0<ct< 3 4 3 d>4, 



(4) 



the BNLS ([2|) is locally well-posed in H 2 . Global existence and scattering of BNLS solutions in 
the i? 2 -critical case a = 4/(d — 4) were studied by Miao, Xu and Zhao }MXZ09j and by Pau- 
sader |Pau09b| . The latter work also showed well-posedness for small data. The _£f 2 -critical defo- 
cusine BNLS was studied by Miao, Xu and Zhao [MXZ08] and by Pausader |Pau071 IP1m09"aj . 

The above studies focused on non-singular solutions. In this work, we study singular solutions 
of the BNLS in H 2 , i.e., solutions that exist in H 2 (fik. d ) over some finite time interval t G [0, T c ), 
but for which lim || "011^2 = oo. The first study of singular BNLS solutions was done by Fibich, 



Ban and Papanicolau [FIP02J, who proved the following results: 

Theorem 1. Let ipo £ H 2 . Then, the solution of the subcritical BNLS ([2|) exists globally in H 2 . 



Then, the solution of the critical BNLS ([3]) exists globally in H 2 . 

The simulations in [FIP02| suggested that there exist singular solutions for ad = 4 and ad > 4, and 
that these singularities are of the blowup type, namely, the solution becomes infinitely localized. 
However, in contradistinction with NLS theory, there is currently no rigorous proof that solutions 
of the BNLS can become singular in either the critical or the supercritical case. 

Most subsequent research of singular BNLS solutions focused on the critical case. Chae, Hong 
and Lee [CHLllJ, showed that radial singular solutions of ([3]) have a power-concentration property. 
In [BFMlOb] . we showed that radial singular solutions are quasi self-similar. We also proved, 
without assuming radial symmetry, that the blowup rate is bound by a quartic-root, the power- 
concentration property and the existence of the ground-state of ©. The two latter properties were 
also proved by Zhu, Zhang and Yang |ZZY10| . In [BFMlObJ, we also provided informal analysis 
and numerical evidence that peak-type singular solutions of the critical BNLS collapse with a quasi 
self-similar R B profile at a blowup rate which is slightly faster than the quartic-root bound. 

In this work, we use asymptotic analysis and numerics to find and characterize peak-type 
singular solutions of the supercritical BNLS. We find that their properties mirror those of the 
supercritical NLS. Ring-type singular solutions of the supercritical BNLS were studied in [BFG10, 



t->T c 




IBFMlOaj . 
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1.1 Summary of results 

We analyze singular solutions of the focusing L 2 -supercritical and -ff 2 -subcritical BNLS, i.e., when 

\A/d<a d<4, 



[A/d<a<^ d>4. 

We assume radial symmetry, i.e., that ip = ip(t,r), where r = |x|. In this case, equation ([2]) reduces 
to 

i^(t,r)-A^ + ^| 2<T V = 0, V(0,r) = Vo(0, (7) 



where 



is the radial biharmonic operator. 

The paper is organized as follows. In Section [21 we show that the supercritical BNLS admits 
the explicit self-similar singular solutions 

where the blowup rate of L(t) is exactly a quartic-root 

L(t) = k{T c - t) 1/4 , k > 0, 
and the self-similar profile S B (p) is a solution of 

-uS B (p) + i^- 05 B + P S B '^j - A 2 p S B + \S B \ 2a S B = 0, 
S B '(0) = 5 B '"(0) = 0, ,S B (oo) = 0. 



(9) 



WKB analysis of the large-/? behavior of shows that it belongs to L 2+2ct , but not to L 2 . Since 
limt-+T c HV'Sb II2+2CT = 00 ' ^ s a singular solution in L 2+2a , but not in H 2 . To the best of our 
knowledge, this is the first time that explicit singular solutions of the BNLS are presented. 

In Section f2. II we show that the zero-Hamiltonian solutions of ([9]) satisfy the boundary condition 

Hm (pS'+(- + i^)s)p'y = 0, Ud-2-1) < 7 <4+^. 
p^oo \ \cjk 4 // 3\ cry a 

In analogy with the supercritical NLS, we conjecture that for any d, a and v, there is unique 
admissible solution S B &dmis ' (p) , which has a zero Hamiltonian and is monotonically decreasing. 
This solution is attained for a unique k = K admis '(a,d,u) > 0. While a rigorous existence proof 
for the S B profile remains open, we provide numerical support for the existence of the admissible 
solutions. 

In Section [3] we consider H 2 singular solutions. Using informal asymptotic analysis and the 
analogy with the supercritical NLS, we conjecture that these solutions undergo a quasi self-similar 
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collapse with the ips B profile, where S B is the unique admissible solution 5 B a mis '. The blowup 
rate of these solutions is given by L(t) ~ K admis (T c — i) 1 / 4 . These characteristics are confirmed 
numerically, in simulations of both the one-dimensional and the two-dimensional BNLS. 

The numerical simulations of the BNLS were performed using the IGR/SGR method [RWOO, 
IDG09| . sec [BFMlOa] for further details. The numerical solution of the nonlinear fourth-order 
ODE for S B is obtained using a modified Petviashvili (SLSR) method, which is described in the 
appendix. The code is available online at http://www.math.tau. ac.il/~fibich/publications. html 

The results of this study are based on asymptotic analysis and numerical simulations, but not on 
rigorous analysis. These results show that there is a striking analogy between collapse of peak type 
solutions in the supercritical NLS and the supercritical BNLS. We note that the rigorous theory 
for singular solutions of the supercritical NLS is much less developed than that for the critical 
NLS. Indeed, a rigorous proof of the blowup rate and blowup profile of the supercritical NLS was 
obtained very recently, and only in the slightly-supercritical regime < ad — 2 <C 1 [MRS09J. We 
hope that this study will motivate a similar rigorous treatment of the supercritical BNLS. 



2 Explicit singular solutions 

Let us look for explicit self-similar solutions of the supercritical BNLS (j2J). Since the BNLS is 
invariant under the dilation symmetry r4 j,i h> > jj/^-ip, where L is a constant, this 

suggests a self-similar solution of the form 

Substituting if)$ B in the BNLS gives 

- r'{t)L\t)S B {p) - iL\t)L t 05 B + psA - A 2 p S B + \S B \ 2a S B = 0. (11) 

Since S B is only a function of p, equation (fTi~|) must be independent of t. Therefore, there exists a 
real constant k such that 

L 3 L t ^ \ (L 4 ) t = -. 4 /4. 

Hence, L(t) is a quartic root, i.e., 



L(t) = Ky/T c -t, k>0. (12) 
Likewise, since 5b is only a function of p, then r'(t)L 4 (f) = v. Hence, 

Substituting (fT2|) and (fl~3|) in ([TT]) shows that the equation for S B is 

- uS B {p) + i^ 05 B + psA - A 2 p S B + \S B \ 2a S B = 0. (14a) 
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Since S B is radially-symmetric and decays at infinity, it should satisfy the boundary conditions 

Sb'(O) = S B "' (0) = 0, 5 B (co) = 0. (14b) 
Equation (|14p has the two parameters v and k. Note, however, that 

S B (p; k, v) := u^Sb^p; k = k/v 1 /*), (15) 
where S B (p; k) is the solution of (|14[) with v = 1, i.e., 



k (2 ~ 



2a _ 



Sn(p) + i'-J [^S B + pS B J- A 2 p S B + S B S B = 0. (16a) 
subject to 

5 B ' (0) = 5b'" (0) = 0, 5 B (co) = 0. (16b) 



Equation (|16p can be viewed as a nonlinear eigenvalue problem with the eigenvalue k and 
eigenfunction S B . By analogy with the supercritical NLS [KL95, BudOl], we make the following 
conjecture: 

Conjecture 3. Let a be in the L 2 supercritical and H 2 -subcritical regime (J6|). Then, there exists 
a solution |5 b (/)),k| to equation (fT6|) . such that S B ^0 and k > 0. 

Hence, we have the following result: 



Lemma 4. Assume that Conjecture holds, and let S B (p; k, v) be a nontrivial solution of ([T] 
Then, 

iPs B (t,r) = t^—S b (-^-] e iuft ^M ds , L(t) = n^/T~~t, (17) 



L 2 /-(t) \L(t) 
is an explicit solution of the BNLS equation (J2j> . 

As p — > oo, the nonlinear term in (|14p becomes negligible, and (|14aj) reduces to 

A 



uS l T{p) - A.pS B n + (~5f> + p(S l B m )p^ = 0, (18) 



2 _ (d-l)(d-3) (d-l)(rf-3) 2 2(d-l) 3 4 



where 



We now use WKB to find the large p behavior of (118p : 
Lemma 5. Let S B n (p) be a solution of (I18p . Then, 

~;lin 



S B n ~ ci5 B)1 (p) + c 2 5 Bi2 (p) + c 3 5 Bi3 (/)) + c 4 5 Bi4 (p), p co, 
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where {cj} 4 =1 are complex constants, 



S B ,i{p) ~ P 

SsM ~ -^^exp(-^M 4 / 3 + z^log(p)), 

^> ~ mp( ;igr 3) ^( + ^ ( - )4/3+ ^H- 

Proof. In order to apply the WKB method, we substitute S^ n (p) = exp(w(p)), and expand 

~ w (p) + wi(p) + ... . 

Substituting u>o(p) = apP and balancing terms shows that p = 4/3, and that the equation for the 
leading-order, the (/o 4//3 ) terms, is 

/ ;\3 ( 4aV V 

H) = y-j) p = l T p ' 

Therefore, 

_3 3 /7^_3k 4 / 3 f . -v/3+i -V3 + il 
The equation for the next order, the (1) terms, is 



, k a (2 A 



implying that 

1 (2 M 4z^ \ , 

«>l = r -(1 - ffd) + — « log p. 
3 \a k 4 / 

The next-order terms are (p- 4 / 3 ) = o(l) and can be neglected. We therefore obtain the three 
solutions jSb 2) S B S , and 5 B 4. 

Since (|18p is a fourth order ODE, another solution is required. To obtain the fourth solution, 
we substitute r^o ~ /31og(p) in ([18]) and obtain that the equation for the leading-order, the 0(1) 
terms, is 

and that the next-order terms are (p -4 ) = o(l) and can be neglected. The fourth solution is 
therefore S B i- 

□ 
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Equation (fl"8j) thus has the two algebraically-decaying solutions, 5 B ,i and S B ,2i the exponentially- 
increasing solution 5 Bj 3, and the exponentially-decreasing solution S B 4. The fact that 5 B) 3 increases 
exponentially as p — > 00 is inconsistent with the boundary condition (|14b|) . Therefore, 

S B (p) ~ ciS B)1 (p) + c 2 S B)2 (p) + C45 , B ,4(p)- (19) 
Since ad > 4, the exponent j^(o~d — 1) of 5 B 2 is larger that the exponent ^ of 5 B) i, hence 

Direct calculations give that 

S B>1 ?L 2 (R d ), 
S Bi2 G L 2 (M d ), 
5 Bj4 E L 2 (R d ), 

Therefore, S B is in L 2+2a . Unless c\ = c 2 = 0, however, S B is not in H 2 . Furthermore, since 

II / ||2+2o- _ 1 || n ||2+2cr 

\\WS B \\2 + 2a - L 4/a-(d-4) U\ H^ 3 H 2+2 CT ' 

then G L 2+2<T for < t < T c . In the F 2 -subcritical regime 4/<r — id — 4) > 0. Therefore, 

t l_™J|V>5 B || 2+2(T = oo. 

Hence, 

Lemma 6. Assume that Conjecture^ holds. Then, ips B is an explicit solution of the BNLS equa- 
tion ([2]) that becomes singular in L 2+2u as t — >• T c . 

2.1 Zero-Hamiltonian solutions 

As is the case of peak- type solutions of the supercritical NLS, a key role is played by the zero- 
Hamiltonian solutions. 

Theorem 7. Let a be in the L 2 supercritical and H 2 -subcritical regime ([6]). Let S B be a solution 
of^M-IfH [S B ] < 00, then H [S B ] = 0. 

Proof. The Hamiltonian of V'Ss , see (|17|) , is equal to 

H \*Sb\ = L 4/.-(,-4) g [Sn] . 

From i? 2 -subcriticality, it follows that L~ A / a ~ i+d {t) 7^ const. Therefore, Hamiltonian conserva- 
tion (H [ips B ] = const) implies that H [S B ] =0. □ 



Sb,1 3* <Sb,2> 



00. 



(20) 



A5 B ,i G L 
AS B<2 L 2 
AS B 4 G L 2 



5 B> i G L 2+2a {R d ), (21a) 
5 Bj2 G L 2+2<T (M d ), (21b) 
S BA G L 2+2<T (IR d ). (21c) 
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Lemma 8. Let a be in the I? supercritical and H 2 -subcritical regime ([6]). Let S B (p) be a zero- 
Hamiltonian solution of (I14p . Then, 



1. c 2 = c 3 = 0. 

2. If ci ^ i/iera 

or, equivalently, 
where 



S B {p) ~ ciS^iO), /)^oo, (22) 



Urn US' + ( \ + ) V ) ,> 0. (2:5a.) 

2 / 2\ 2 

7o<7<7l, 7o=o d ~ 2 -- > 71 = 4 + -. (23b) 

3 V <t / a 



Moreover, S B G L 2(J + 2 and S B £ L 2 . 

Proof. The exponentially increasing solution 036^,3 must vanish, as explained above. Convergence 
of the Hamiltonian requires that AS B G £ 2 . Since AS B 2 $ L 2 , see f)21 [) . it follows that C2 = 0. 
bmce S B 1 L 2 , then Sb L 2 . 

To show that (|23|) is equivalent to demanding that C2 = 0, we first note that in the ff 2 -subcritical 
regime d — 4 < - and so 70 < § (2 + ^) < 71. Next, direct calculation gives that 



00, 



and that 

where the LHS is the result of the next term in the WKB approximation of S B i- Therefore, 



P> 



P T + I + ^ ) Sb ~ ° ( Cl ' P "~°~ 4 ) + °( C2 ' ^ 7+4/3 -^ M - 1} ) , P^oo. 



dp a k 
Since, for 70 < 7 < 71, 

lim pT-l" 4 = 0, lim = 00, 

p— >oo p— >oo 

it follows that C2 = if and only if the limit (j23aH is satisfied. □ 

The fourth-order nonlinear ODE (|14ap requires four boundary conditions. Three boundary con- 
ditions are given by (|14bj) . and the fourth condition will be the zero-Hamiltonian condition (|23ap . 
Generically, one can expect that for a given u, this nonlinear eigenvalue problem has an enumer- 
able number of eigenvalues with corresponding eigenfunctions S B ^ n \ As in the case of the 
supercritical NLS [KL951 iBudOl] . we conjecture that for any (er, d, v) there is a unique admissible 
solution, which is monotonically decreasing. 
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Conjecture 9. Let a be in the L 2 supercritical and H 2 -subcritical regime and let v > 0. Then, 
the nonlinear eigenvalue problem posed by equation (j!4a|) . subject to the boundary conditions 

S B '(0) = SJ"(0) = S B (oo) = 0, lim f pS' +(- + iK] s) p* = 0, (24) 
where 7 satisfies (|23bj) . admits a unique eigenpair (5 B admis (/9),K admis -), such that 

,.admis. admis./ j , \ ^ n 

k = k [a, a, v) > U, 
and I ,S' fl admis - (p) J i s monotonically- decreasing. Furthermore, 

q admis. / \ „—2/a—4ui/n 4 „ . 
S B (p) ~ cp 7 7 , p ► OO, 

and 

.^admis./^ j ,,\ , ,1/4 -admis. /_ j\ -admis. . ,^admis./„ j i \ 

k (<r, ct, z^J = v 1 k [a, a), k := k {a,a,i' = 1). 

3 Peak-type i7 2 -singular solutions 
3.1 Informal analysis 

As in the supercritical NLS [LPSS881 ISKL93| . we expect that singular peak-type solutions of the 
supercritical BNLS undergo a quasi self-similar collapse, so that 

W,r)~(^ (t ' r) 0<r<r c , ^ 

I ^non-singular (^> T ~>_T C , 

where ips B * s the self-similar profile (|10p . The singular region r £ [0, r c ] is constant in the coordi- 
nate r. Therefore, in the rescaled variable p = r/L(t), the singular region p € [0, r c /L(t)] becomes 
infinite as L(t) — > 0. This is in contradistinction with the critical-BNLS case, where the singu- 
lar region p G [0, p c ] is constant in the rescaled variable p, but shrinks to a point in the original 
coordinate r |BFM10b] . 

Lemma 10. Let ad > 4, and let tp be a peak-type singular solution of the BNLS that collapses with 
the i^Sb profile (fTUj) . If L(t) ~ k(T c — t) p , then p>\. Furthermore, 

• If P = 1/4 if then the self- similar profile S B (p) satisfies the equation (|14a|) . 

• If p> l/A, then the profile satisfies the equation 

-A 2 R B (x.)-R B + \R B \ 2r7 R B = 0. (26) 

Proof. If ij) ~ "05b, then the equation for S B is 



uS B - i (lim L t L^j + P S B '^j - A 2 p S B + \S B \ 2a S B = 0, 



(27) 



implying that LtL 3 should be bounded as t — > T c . Since L 3 Lt ~ —pK^{T c — t) Ap 1 , it follows 
that p > 3. If p = 1/4, then equation ([27|) reduces to (|14aj) . see Section [2j □ 
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From Hamiltonian conservation it follows that H [ips B ] 1S bounded, because otherwise the non- 
singular region would also have an infinite Hamiltonian. Therefore, from Theorem [7] it follows 
that H [S B ] = 0. 

In Lemma [8] we saw that the zero-Hamiltonian solutions of (|14ap are in L 2+2a , but not in L 2 . 
Hence, ips B L 2 . From power conservation, however, it follows that if ipo E H 2 , then ^ G L 2 . As 
in the NLS case, sec [BP92], this "contradiction" can be resolved as follows. 



Corollary 11. Let S B (p) be a zero-Hamiltonian solution of (|14p . Then, \\S B \\ 2 = oo. Neverthe- 
less, Hm Us B \\ L 2 {r<rc) < oo. 

Proof. Since 5 B (p) ~ ci5 B i(p), 



/■oo 

ll^illa ~ c / P" 4/<7+d_1 dp ~ C/~ 4/<J 

Jp=0 

The profile V'S'b satisfies 



oo 

= OO. 

p=0 



r r c/L(t) 

mW< m = ^- 4/ct (o- / isbWiv- 1 ^ 



p=0 



L d ~^(t) ■ {C(M^^\ = 0(1). 



□ 



In summary, we conjecture the following: 
Conjecture 12. 

Lei V ^ e peak-type singular solution of the supercritical BNLS. Then, 

1. The collapsing core approaches the self- similar profile ips B > *- e -> 

1>{t,r) ~ ilJSB&r), 0<r<r c , (28a) 

where 

2. The self-similar profile S B (p) = S B admis ' (p) is the unique admissible solution of 

.k* (2 



vS B (p) + ( -S B + P SJ - A 2 S B + \S B \ 2 °S B = 0, 

4 \a J H (28c) 

S B '(0) = S B '"(0) = 0, S B (oo) = 0, H[S B ] = 0, 

where K admis -(a,d,v) = v^^^ia.d) > 0. 
3. In particular, S B (p) ^ R B {p)- 
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\\\r\ 



0.001 0.002 
t 




(A) d = 1, a = 6, Vo = l.6e~" (B) d = 2, a = 3, Vo = 3e" r 

Figure 1: Maximal amplitude of singular solutions of the supercritical BNLS. 



4- The blowup rate of singular peak-type solutions is exactly a quartic root, i.e. 



L(t) ~ Ky/T c -t, k>0. (28d) 

5. The coefficient k of the blowup rate of L(t) is equal to the value of k of the admissible solu- 
tion S B , i.e., 

k := lim = K admis -(cj,(i,i/). 

In particular, k is universal (i.e., it does not depend on the initial condition). 
In Section HOI we provide numerical evidence in support of Conjecture [T2l 

3.2 Simulations 

The radially-symmetric BNLS ([7]) was solved in the supercritical cases: 

1. d= 1, a = 6 with the initial condition ipo( x ) = l-6e _:r . 

2. d = 2, a = 3 with the initial condition ijjo{r) = 3e _r . 

In both cases, the solutions blowup at a finite time, see Figure [TJ 

To check whether the solutions collapse with the self-similar profile (|28p. the solution was 
rescaled according to 

VWaied (t, P ) = L 2 ^(t)^(t,r = p-L), L(t) = ||V||~ 2At . (29) 

Comparing this rescaling with (|28b|) shows that it implies that 1^(0)1 = H-SbH^ = 1- This require- 
ment can always be satisfied with a proper choice of u, see (|15p . Figure l2Al shows the rescaled 
solutions at the focusing levels L = 10~ 4 and L = 10 -8 , and the rescaled solution S B (p) of (|28|) F*1 



1 In the calculation of 5b, see Appendix [XJ the values of v and k were extracted from the BNLS simulation as 
discussed below, see equations (|30I32[) . 
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L \\\t\ 



r 2/o i i 



'^rl_ — j 1 nl i 1 

5 10 15 5 10 15 

x/L r/L 

(A) d= l,cr = 6 (B)d = 2,cr = 3 

Figure 2: The solutions of Figured! rescaled according to (|29p . at the focusing levels 1/L = 10 
(blue solid line) and 1/L = 10 s (black dashed line). The red dash-dotted line is the rescaled 
solution 5 B of (f28|) . The magenta dotted line is the rescaled ground-state R. 




(A) d= l,cr = 6 (B)d = 2,(T = 3 

Figure 3: The solutions of Figure [H rescaled according to (f29|) . at the focusing level 1/L = 10' 
(circles). Solid lines are the fitted curves y = 0.63 • (x/L)-°- 33 (left) and y = 0.85 • (r/L)-° M< 
(right). 
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(A) d= l,cr = 6 (B)d = 2,cr = 3 

Figure 4: Convergence to a self-similar profile. The solutions of Figure HJ rescaled according 
to ([29]) . as a function of log(r/L), at the focusing levels L = 10 _1 (dashed blue line), L = 10 -2 
(dash-doted red line), L = 10~ 3 (dotted green line), L = 10 -4 (solid black line) and L = 10 -8 
(solid magenta line). The circles mark the approximate position where each curve bifurcates 
from the limiting profile, see also Table [TJ 



l/L 


10 


100 


1000 


10000 | 


r c 


x/L (d = 1) 


3.6 


36 


360 


3600 


0.36 


r/L (d = 2) 


6 


60 


600 


6000 


0.6 



Table 1: 



'osition of circles in Figure El 



The three curves are indistinguishable, showing that the solution is self-similar with the S B profile, 
and not with the R B profile. As additional evidence, Figure [3] shows that as p — > oo, the self-similar 
profile of ip decays as p~ 2 ^ a , which is in agreement with the decay rate of S B i(p). 

We next verify that the solution converges to the asymptotic profile for r £ [0, r c ], i.e., for p £ 
[0,r c /L(t)\. To do this, we plot in Figure H] the rescaled solution at focusing levels of l/L = 
10,100,1000,10000, as a function of log(r/L). The curves are indistinguishable at r/L = 0(1), 
but bifurcate at increasing values of r/L. These "bifurcations positions" are marked by circles in 
Figure HI and their r/L values are listed in Table [TJ The "bifurcation positions" are linear in l/L, 
indicating that the region where ip ~ Vs B ^ s indeed p E [0, r c /L(t)], which corresponds to r £ [0, r c ]. 

In order to compute the blowup rate p, we performed a least-squares fit of log(L) with log(T c — t), 
see Figure [5j The resulting values are p 0.2502 in the d = 1, a = 6 case, and p « 0.2504 in 
the d = 2, a = 3 case. Next, we provide two indications that the blowup rate is exactly 1/4, i.e., 
that 

L(t) ~ k{/T c - t, k>0. 
1. If the blowup rate is exactly a quartic root, then L 3 Lt — > < 0. Indeed, Figure l6Al shows 
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(A) d= l,a = 6 
Figure 5: L(t) function of (T c 

of Figure [Q (circles). Solid lines are 
and L = 0.931 • (T c — i) - 2504 (B). 




(B) d = 2, a = 3 
- i) , on a logarithmic scale, for the solutions 
the fitted curves L = 1.048 • (T c — t) - 2502 (A) 




(A) 



0.4- 



lV 



0.2 



10 10 



10' 

1/L 



10 10 



(B) 



Figure 6: A: L 3 L^ as a function of 1/L, for the solution of Figure [TAl (black solid line) and of 
Figure [TBI (red dashed line). B: Same as (A) for L 4 Tt, where r = argip(t,r = 0). 
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that in the case d = 1, a = 6, L 3 L^ — > —0.289, implying that 

K admis. (d = i, = 6 ) « ^4 . 0.289 ss 1.037 . (30a) 

In the case d = 2, <r = 3, L 3 Lt — > —0.171, implying that 

K admis.( d = 2) ff = 3) ra ^ 4 . 0-171 ~ 0909 _ ( 30b ) 

Since L 3 Lj converges to a finite, negative constant, this shows that the blowup rate is ex- 
actly 1/4. 

2. According to Lemma [TUl if lim.t-+T c L Lt < 0, the self-similar profile S B (p) does not satisfy 
the standing- wave equation ([26]) . but rather is a solution to the problem (|28c|) . as is clearly 
demonstrated in Figure [2j 

We recall that calculation of the profile S B requires knowing the numerical values for v and k, 
see Appendix |Aj The value of n was obtained previously from the limit limi_>.T c (— 4L 3 Lj) 1//4 . We 
approximate the value of the coefficient v from 

dr 

^numeric = lim L 4 (t) — , t = argx(j(t,r = 0). (31) 
t^T c dt 

Indeed, Figure l6Bl shows that in both cases L 4 (t)^- quickly converge to 

^numoric(d = 1, a = 6) = 0.36187, z. numeri c(d = 2, a = 3) = 0.22826. (32) 

As a further verification, using the above values for v, we seek a value of k such that the solution 
of (I28cp will satisfy |S' B (0)| = 1, and obtain 

K (d = 1, a = 6, v = 0.36187) = 1.007, n{d = 2, a = 3, v = 0.22826) = 0.894. 

These values are within l%-3% from the values of 1.037 and 0.909 we obtained directly from the 
BNLS simulations, see ([30]) . 

Finally, we verified that the value of k in the blowup rate (|28dj) is universal. We solve the BNLS 
in the case d = 1, a = 6 with the initial condition ipo(x) = 2e~ x . In this case, the calculated value 
of n(d = l,a = 6) is k = lim v^— 4LjL 3 ~ 1.037, which is equal, to first 3 significant digits, to 

t^T c 

the previously obtained value, see (|30ap . for the initial condition ipo(x) = 1.6e~ x . Similarly, in the 
case d = 2, a = 3, we solve the equation with the initial condition ipo(x) = 3e~ x . The calculated 
value of n{d = 2, a = 3) is k = lim y— 4L^L 3 ~ 0.913, which is equal, to first 2 significant digits, 

t^T c 

to the previously obtained value, see (|30bj) . for the initial condition ipo(x) = 3e~ x . 
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A Numerical calculation of the S B profile 



In order to solve equation (|28cp . we first define its linear part, which is the fourth-order linear 
differential operator L [S] 



L[S(p)] = -vS{p)+i^-(^S + pS^-tf p S, 
under the BCs, see equation 



(33a) 



S'(0) = 5"'(0) = 5(oo) = 0, lim pi ( pS' + (- + A ) S 



70 < 7 < 7i; 7o 
The nonlinear ODE ()28c|) is therefore rewritten as 



2 / 2 
-U-2-- 

3 V cr 



(33b) 



7i = 4 + -. 

a 



L[S{p)] + \S\ 2a S = 0. 



(34) 



For given numerical values of v and k, we wish to calculate the ground state of the nonlinear 
boundary- value problem ()34p . In order to do so, we modify the SLSR method for the calculation 
of the ground-state of the NLS |Pet76l [PSOT] IAM05j and BNLS [BFMlObj as follows. We consider 
the fixed-point iterative scheme 



(p) 



-L 



-i 



k = 0,1,. 



(35) 



for the solution of (I34p . In the standard application of the SLSR method, L is a differential operator 
of constant coefficients, and its inversion is easily performed using the Fourier transform. In our 
case, L is a variable-coefficient operator, and the Fourier Transform cannot be used. Therefore, 
we discretize the operator L using finite differences, see Appendix IA.1} and invert it using the LU 
decomposition. 

We observe numerically that generically, the iterations (|35p converge to zero for a small initial 
guess and diverge to infinity for a large initial guess. To avoid this divergence, we rescale the 
approximate solutions at each iteration, so that they satisfy the integral relation: 



\S\ 2 p d - 1 dp= (S,S) = —Re(S,L~ 1 \Sf Sa S) 



ll al2o- , 



which follows from multiplication of (J35J) by S. Here, < • > denotes the standard inner product 
(/><?) = / f* 9P d l dp. Following a similar argumentation as in [BFMlOb] . we obtain that the 
iterations are 



.(5« )(S (fe)) 



Re(S( k ),L- l \S( k )\ 2 °S( k )) 



(36) 
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In our simulation, this method converged for every value of v and k that we tried. The numerical 
values of v was obtained from the on-axis phase of the BNLS simulation solutions, as explained 
in Section 13.21 In order to obtain a prediction of k, we recall that the specific choice (I29D of the 
blowup rate L(t) implies that |5 B (0)| = 1. We therefore use the SLSR solver to search for the value 
of k for which | S B (0) = 1 1 . 



A.l Discretization of L 

Using half-integer grid 



Pn 



n + - ] h, 



n 



0,...N-1, 



R, t 



N 



the O (/i 2 ) centered-difference discretizations of the radial biharmonic operator D 2 p and of the first- 
derivative, the approximation at the interior nodes is 



IK 



IK 



~>n+l 



?n-l 



2h 



DpS n + \S n \ S n 



Oih 2 



The stencil is five-nodes wide, so two ghost-nodes are needed at each boundary. In order to 
enfold the ghost-nodes at p = 0, we relate them to the interior nodes, using the symmetry of the 
solution S(p) = S(—p), so that 



" S_ 2 " 




" 1 " 




" So ' 


. S-l _ 




1 




. Si _ 



This relation is substituted in the discretization of the equation at po and p\ . 

At the other boundary p = i? max we use the approximate form of the solution obtained from 
the WKB approximation, i.e., we require 



S n = ClS B<1 (p n ) + C 4 S , b i4 (/)„), 

In matrix form, this becomes 

S b ,i(pn-i) Sn t 4{pN-i) 

Sb,i(pn) S b ^(pn) 

Sb,i(pn+i) S B)4: (pn+i) 

Sb,i(pn+2) S B)4: (pn+2) 

which is then solved to obtain 



Sb,i(pn+i) S B)4: (pn+i) 
Sb,i(pn+2) S B)4: (pn+2) 



n 



N - 2,N - 1,N,N + 1, 



Sn+i 




Sn+2 











Sn-i 




Cl 




Sn 




C 4 _ 




Sn+i 








Sn+2 


B,l(PN- 


-l) 


SB,4:(PN- 



Sb.i(pn) S b a(pn) 



-1 


Sn-i 




Sn 



Some care should be taken when choosing the parameters -R max and N. On the one hand, we 
use the closed- form approximations for S B i and S B 4 that become more accurate for -R max 3> 1 • 



REFERENCES 



18 



On the other hand, since S B ,4 has a super-exponentially decreasing term e p , choosing too- 
large a value of -R max leads to numerical instabilities. Finally, in order to resolve the rapid- 
oscillations e ip4/A of 5 B ,4, the grid-size N must be chosen such that p l ^h p = O (^Rmax/N^ <C 1, 

hence that N > O ^-Rmax) • The grid-size N , however, cannot be arbitrarily large, since the con- 
dition number of L is O (-ZV 4 ) . 

In the simulations presented in this study, we used an extension of above approach to a fourth- 
order approximation, and set R max = 160 and N = 32000. 
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